Graphene nanodevices: bridging nanoelectronics and subwavelength optics 
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The unconventional properties of graphene, with a massless Dirac band dispersion and large 
coherence properties, have raised a large interest for applications in nanoelectronics. In this work, 
we emphasize that graphene two dimensional character combined with current standard lithography 
processes allow to achieve devices smaller than the Dirac electrons wavelength. In this regime, 
we demonstrate that the electronic properties present deep analogies with subwavelength optics 
phenomena. We describe the rich transport physics in graphene-based nanodevices through optical 
analogies: From the Bethe and Kirchhoff-like diffraction patterns in the conductance of graphene 
slits to the Fabry- Perot oscillations of the conductance in nanoribbons. We introduce the concept 
of electronic diffraction harriers^ which transmission cancels at the Dirac point. This gives central 
insight in the properties of Graphene subwavelength devices including nanoelectronics standard 
systems, such as quantum dots. As an application we propose a new type of quantum dots, namely 
functionalized subwavelength quantum dots, which could be used as molecular spin valves. 



Analogies play a prominent role in physics. By allow- 
ing transfer of notions and advances from one field to 
another, they can provide deeper insight in those fields. 
In this context, the analogy between optics and quantum 
electronic transport, thanks to the ondulatory nature of 
electrons and light, has been particularly useful either in 
predicting the transport properties of devices [1 , or in 
understanding optical phenomena such as the coherent 
multiple scattering of light [2j. Indeed, as long as the co- 
herence lengths of a material are large, the -quantum 
transport- Landauer formula establishes a one to one 
correspondence between conductance and transmission 
in scattering experiments. While light scattering exper- 
iments probe the transmission probabilities between a 
given pair of input and output modes, the transport ex- 
periments probe the sum of the transmission probabili- 
ties from all possible input modes to all possible output 
modes at the Fermi energy, i.e. at a given frequency in 
the optical language. 

Such analogies are particularly relevant in graphene - 
and have already been exploited [3]- due to extraordinary 
coherence properties of this material. Indeed, this system 
offers unique characteristics either in the exfoliated or in 
the epitaxial form [4[ |5] IH [J • The electronic mean- free 
path can be of the order of a micron at room temperature 
and the Dirac electrons wavelength, which in principle 
diverges at the Dirac point, can be up to 100 nm or even 
more in real systems. Moreover, graphene is at the 
surface and directly accessible to lithography, allowing 
to produce nanodevices down to a few nanometers as 
typical size [9l [10], i.e. smaller than the Dirac electron 
wavelength. 

In this work, we demonstrate that the transport prop- 
erties of such devices can be understood by bridging na- 
noelectronics of graphene nanodevices with the domain of 
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subwavelength optics. Indeed, at this scale, quantum in- 
terference effects in graphene become prominent, and the 
conductance exhibits a universal behaviour characteris- 
tic of subwavelength optics. We believe that this analogy 
between two fields, where intense effort are made to find 
devices with original properties, can be very fruitful and 
we give examples in this work. 

We consider first systems such as slits and nanoribbons 
sandwiched in between graphene semi-infinite sheets as 
leads, which have well known optical analogues. These 
studies lead us to the mixed "electronics-optics" concept 
of electronic diffraction harriers which transmission can- 
cels for a divergent wavelength of the electron in the 
graphene plane, i.e at the Dirac point. This concept 
gives also an important insight in the confinement in 
graphene quantum dots where a chaotic Dirac billiard 
behaviour has been recently observed [11 . It shows that 
the characteristics of these quantum dots, such as barrier 
transparencies, can be tuned and controlled, thus open- 
ing to a wealth of new potential applications. A new type 
of quantum dot is proposed by functionalizing carbon 
atoms in a graphene constriction. In particular function- 
alization by magnetic molecules could provide molecular 
spin valves. 

The prototypical system we consider is shown in Fig.[T] 
It is a tripartitioned system consisting of a central device 
(a graphene nanoribbon, a slit or a graphene quantum 
dot), contacted to mesoscopic leads made of two semi- 
infinite graphene half planes. The system parameters 
which affect the response of the system with respect to 
optical-like phenomena are the aperture W of the con- 
tact, the typical length L, and the energy dependent elec- 
tron wavelength A. The chirality and the electronic struc- 
ture of the ribbon can also affect the conductance of the 
device. Yet we will show that the quantum interference 
effects studied here present a universal behaviour, inde- 
pendent of the electronic structure of the system, that 
present a deep analogy with phenomena that are known 
in subwavelength optics [T2j [13] . 
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FIG. 1: Schematics of the systems studied in this article. Left: 
the yellow part represents the metallic contacts and the blue 
region is the graphene part of the device. Right: detailed 
view of the tripartitioned system. The typical width of the 
contact is W and the typical length of the central device is L. 
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FIG. 2: Conductance of graphene slits and diffraction Left: 
schematic view of the geometry of a graphene slit. Right: 
Scaling law for a graphene slit conductance as a function of 
W/X. This shows the Bethe-like X ^ W quadratic regime 
and the Kirchhoff-like X <^ W linear regime (see text). The 
inset shows the conductance of slits with different widths. 



Our methodology relies on an exact numerical calcula- 
tion of the conductance within a tight-binding and a 
Landauer quantum transport formalisms. The novelty is 
that the Landauer equation is solved via a new recursive 
numerical algorithm which reduces complex 3D/2D var- 
ious shaped quantum transport devices into an effective 
ID system. Here only the channels effectively contribut- 
ing to the conductance are considered; the method thus 
allows to exactly compute the contact resistance in an 
efficient way. 

Conductance of graphene slits and diffraction - The 
archetype experiment showing an optical behaviour is the 
transmission through a hole, for example in the geometry 
indicated in the left panel of Fig. [2] The conductance of 
the slit is related to the transmission of eigenstates (or 
modes) through the slit. Therefoore the situation is obvi- 
ously connected to the textbook experiment of diffraction 
by a hole. In optics depending on the hole width relative 
to the wavelength there are different regimes. The stan- 
dard Kirchoff theory applies only in the large wavelength 



limit. However the phenomenon of diffraction through a 
subwavelength slit or through an array of such slits is cen- 
tral to subwavelength optics and plasmonics [121 [131 [E] • 
Analogeous results are detailed here forthe conductance 
of graphene slits. 

Here we report the conductance of several slits that 
consist of a single motifs i.e. a single hexagon (the edges 
of the graphene half-planes being of the armchair type) 
but differing by the width W (Fig. [2] inset). We introduce 
the wavelength A of the incident Dirac electrons: 

271 h hvp , . 

where k is the Dirac electrons wavevector, p is the mo- 
mentum, E = vyP is the energy, linearly dispersed in 
graphene with a Fermi velocity oi — 10^ ms~^ [16j. 

We first note (Fig. [2] inset) that for all the slit widths 
W a universal scaling occurs such that the conductance 
depends only on the ratio W/X. This scaling breaks down 
when the wavelength is A ^ 4 nm, that is the energy is 
greater than approximately 1 eV. Indeed as shown below 
this scaling is intimately related to the scaling properties 
of the Dirac equation which is valid only in the low energy 
or long wavelength limit. 

We then note that in the limit W/X <C 1, the con- 
ductance is quadratic, with the conductance g{W/X) 
25(1^/A)^, while in the opposite limit W/X > 1, the con- 
ductance turns out to be linear with g{W/ A) :^ 2.b{W/ A) . 
The crossover between the two regimes occurs around 
W/X :^ 0.3. This promptly reminds of the analogy with 
classical optics and hence offers an immediate interpre- 
tation. We can identify two different diffraction regimes: 
for wavelengths X ^ W^ much larger than the aper- 
ture, we observe a Bethe-like diffraction regime [17] with 
a slit transmitted energy proportional to the square of 
the aperture. On the other hand, for A <C one has a 
Kirchhoff-like diffraction regime^ the transmitted energy 
being proportional to the aperture. Thus the fact that 
the system does not conduct at the Dirac energy is a 
manifestation not of the system electronic structure, but 
rather of the classical optical inability of an incident wave 
to be transmitted through a subwavelength constriction 
which hence behaves like a barrier. 

The explanation of both the observed diffraction pat- 
terns and of the universal scaling law comes from analo- 
gies between Maxwell equations and the Dirac equation 
here. Indeed, a universal scaling is also observed in ordi- 
nary optical diffraction for the slit transmitted intensity 
due to the scaling invariance of Maxwell equations. Here 
it stems from the fundamental scaling invariance of the 
graphene 2D Dirac equation which can be written in a 
dimensionless form, 

nv^a'^=E^ =^ t7.|^=2^V^, (2) 

where a are the Pauli matrices, is a two component 
pseudo-spinor, and f = r/A. The scale invariant (r/A de- 
pendent) form of the Dirac equation directly stems from 
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the linear dispersion relation E = hv^/X. This leads to 
the fact that the conductance undergoes a scale invari- 
ance g{xX^xW) = g{X^W)^ so it only depends from a 
reduced argument W/X^ g{X^W) = g{W/X). 

The argument is as follows. Let us consider a sys- 
tem where each half plane of graphene is replaced by 
a ribbon of width W. The symmetry axis of the rib- 
bon is perpendicular to the slit, passing by the center 
of the slit. The conductance for half graphene plane is 
g{X,W) — lim^^^G{X,W,W). We assume that the 
boundary conditions preserve the scaling, which is evi- 
denced here by the numerical results. This is also the case 
for continuous type models of confinement [T8^. Then 
there is a one to one correspondence between the scat- 
tering states at wavelength A ribbon width W, slit width 
W and those at wavelength xA, ribbon width xW and 
slit width xW. This means that the conductance satis- 
fies G(A, W) = G(xA, xW^ xW) and as a consequence 
g{xX,xW)=g{X,W)=g{W/X). 

Nanoribbons as subwavelength waveguides - Another 
very well-known evidence for interference effects, are the 
Fabry-Perot oscillations, which occur when an optical 
waveguide is sandwiched in between two reflecting sur- 
faces. Several studies of subwavelength optical waveg- 
uides have shown the existence of Fabry-Perot oscilla- 
tions of the transmitted intensity for an incident wave 
parallel to the axis of waveguide. This is for instance the 
case for strips in subwavelength optical metallic gratings 
[19] or subwavelength slits in metallic screens with a fi- 
nite thickness [20] [21] . In these systems the refiection at 
the ends of the waveguide is in fact intimately related to 
the diffraction by the subwavelength aperture. 

In our case, the corresponding 2D transport config- 
uration is a metallic entity (corresponding to the opti- 
cal waveguide), that we chose to be a metallic graphene 
nanoribbon, connected to two semi-infinite planes of 
graphene. Here also as discussed below the diffraction by 
the small apertures of the ribbon produces the refiection 
phenomena at the ends of the nanoribbons. Graphene 
nanoribbons can be realized by standard lithography 
techniques [9l [TOl [23j. We considered a zigzag and an 
armchair ribbon, both chosen with a metallic character 
and hence an available conductance channel at the Dirac 
energy [22 , which can thus be considered from an optics 
point of view as waveguides. 

Calculating the conductance for both metallic nanorib- 
bons (Fig.|3| clearly reveals a Fabry-Perot behaviour. Let 
us note that these oscillations exist for perfect nanorib- 
bons but could be destroyed for example by edge defects 
[23( l24( |2^ . The most evident features are large ampli- 
tude oscillations in the conductance, from maxima val- 
ues of 1 to minima placed along an envelope function. 
These are standard Fabry-Perot oscillations of the trans- 
mittance which follow the Airy function 
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FIG. 3: Graphene nanoribbons as subwavelength waveguides 
Left: schematic view of nanoribbons contacted with graphene 
half-planes. Right: Fabry-Perot oscillations of the conduc- 
tance of armchair nanoribbon with a length of 6 nm (upper 
part) and zigzag nanoribbon with a length of 3 nm (lower 
part). The envelope of the minima of the conductance os- 
cillations deduced from Eq. Q is represented by the dashed 
line. 



where (j){E) = 2k{E)L + 2^{E) is the phase difference 
after one loop in the nanoribbon. k{E) is the wavevector 
of the Bloch state in the infinite ribbon with energy E. 
(j){E) is the phase factor acquired at each refiection. L 
is the length of the Fabry-Perot interferometer. F(yE) = 
m{E)/{l - R{E)f is the finesse coefficient where R{E) 
is the reflection coeflicient at each end of the nanoribbon. 

The Airy function presents maxima = 1 equal 

to 1 when 0/2 is an integer m multiple of tt. For sufli- 
ciently large L the phase (t){E) = 2k{E)L + 2^{E) varies 
rapidly with energy, as compared to F{E); and the min- 
ima occur essentially when sin^((/)/2) is maximum, that 
is at (j)/2 = mir + 7r/2, and envelope along the function 
T^^{E) which independent of the size as we have nu- 
merically checked. 



Timin / 77i\ 
FP [^) 



1 



^ {i-R{E)r 

1 + F{E) (l + i?(^))2- 



(4) 



TME) 



1 + F{E) sin%(p{E)/2) 



(3) 



The minimum of conductance tends to zero at small en- 
ergies indicating that R{E) tends to 1 close to the Dirac 
energy. 

The peaks full width at half maximum is ^^fwhm = 
2{1 — R)/ \/R such as the peaks look thin when R is close 
to 1, and broaden when ^ 0. In our case when E ^ {) 
the peaks seem to have negligible width, indicating again 
that R ^ 1 close to the Dirac energy. The peaks broaden 
for the higher energies, indicating lower values of R. 

Concept of electronic diffraction barrier - The previous 
study on nanoribbons shows the existence of Fabry-Perot 
oscillations and thus a phenomenon of reflexion at each 
end of the ribbon characterized by a reflexion coefficient 
R{E). Just as for the optical analogues this relexion 
phenomena is intimately connected to diffraction at the 
ends of the ribbon. 

The reffexion coefficient R{E) is related to the con- 
ductance G of a junction between a semi-inflnite metal- 



4 




Energy [eV] 



FIG. 4: Electronic diffraction barrier at a contact. Left: the 
contact geometries of the armchair (upper panel) and zigzag 
(lower panel) semi-infinite ribbons coupled to graphene half- 
planes. Right: Conductance of the device represented on the 
corresponding left part, g = T2e^ /h with T the transmittance 
of the diffraction barrier. 

lie ribbon and a graphene half-plane (Fig. [4| by G = 
{2e^ /h) X (1 — R[E)). If there was no reflexion we would 
expect to observe a conductance equal to 1 quantum of 
conductance (2e^ /h) (dot-dashed line in Fig. [4| in the 
range of the Dirac/Fermi energy. Instead, the calculated 
transmittance T{E) = 1 — R{E) through the junction 
(continous line in Fig. [4| is observed to drop to zero at 
the Dirac energy, no matter the chirality, exactly like in 
slits. This means that the reflexion coefficient R{E) tends 
to 1 at the Dirac energy as expected from the analysis 
of Fabry- Perot oscillations (see above). We have checked 
also that formula |4] gives the correct minima of the os- 
cillations with R{E) calculated from the junction trans- 
mission T{E) = 1 - R{E). 

Let us stress that T 7^ 1 behaviour is genuinely a man- 
ifestation of the contact resistance at the junction, which 
is independent of the length of the ribbon, due to its 
metallic character. At the opposite, in the case of non- 
metallic nanoribbons where confinement effects create a 
gap in electronic structure P [7] , the average resistance 
increases exponentially with the length of the ribbon. 

The junction conductance we have calculated here is an 
example of the elementary electronic diffraction barrier 
which is at the basis of the universal behaviour observed 
in the devices hereby studied. 

Application to quantum dots - With the concept of elec- 
tronic diffraction barrier at a contact we can now analyse 
the conductance of standard systems of nanoelectronics 
i.e. quantum dots. We consider two examples that il- 
lustrate the role of geometry and the role of a chemical 
functionalisation of the dot . 

We consider first a graphene quantum dot (Fig.[5| con- 
sisting of a purposely irregular shape graphene nanos- 
tructure contacted via small apertures to the two half- 
planes graphene leads. From the concept of electronic 
diffraction barrier we expect that the quantum dot is 
weakly coupled to the leads when the energy of the elec- 
trons in the graphene sheet is close to the Dirac energy, 
i.e. when their wavelength is sufficiently large. 




FIG. 5: Effect of the geometry on the conductance of quantum 
dots showing well defined resonances close to the Dirac energy. 
Inset: geometry of the irregular shaped quantum dot coupled 
to the graphene half-planes. 




FIG. 6: Effect of a chemical functionalization on the conduc- 
tance of quantum dots Left: geometric configuration with the 
different post ions of the functionnalized atom. Right: con- 
ductance curves corresponding to the different positions of 
the functionnalized atom. The blue curve is the reference 
conductance i.e. before functionalization. 



We note that the conductance of this system (Fig. [5| 
presents a maximum around 3 eV which is the energy of 
the Van-Hove singularity of bulk graphene where the den- 
sity of states diverges. We also note many irregular peaks 
superimposed to the envelope, clearly showing character- 
istic resonances of the quantum dot to be associated to 
the particular (irregular) shape. The conductance de- 
creases when approaching to the Dirac point, which is a 
signature of the lowest number of available states and of 
a diffraction barrier at the contact constrictions. 

The maximum to minimum transmission ratio between 
adjacent peaks is much larger in the vicinity of the Dirac 
energy than at other energies, a phenomena which is also 
observed in the nanoribbons studied previously. This is 
due to the two electronic diffraction barriers at the con- 
tacts which are nearly completely reflecting {R = 1) at 
the Dirac energy. This leads to well defined states (or 
modes) within the quantum dot, and thus to transmission 
channels at well defined energies through these states. 

A second example is the effect of the functionalization 
by impurities on the response of a graphene quantum dot. 
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Here our quantum dot is a very short graphene nanorib- 
bon contacted by small apertures to half-plane graphene 
leads. This is like in the slit geometry studied above. The 
novelty is that now resonances in the quantum dot are in- 
duced by functionalization. A carbon atom is removed at 
a given position and this should mimic the effect of func- 
tionalisation by chemical entities like OH or other groups 
that strongly hybridize with the carbon Pz orbitals [26 . 
Due to the specific electronic structure of graphene, func- 
tionalization of one carbon atom in the plane creates a 
so-called vacancy state at the zero {i.e. Dirac) energy 
[27] . This vacancy state deeply modifies the conduction 
properties of the constriction. Indeed Fig. [6] shows at the 
Dirac energy a well defined peak of conductance corre- 
sponding to the vacancy state. 

When the defect is applied within the short nanorib- 
bon, like in the positions indicated as a), c) or d) in 
Fig. [6} the transmission peak reaches the full value 1. 
The total conductance characteristics is also affected by 
the lateral position of the defect. A defect whose posi- 
tion is shifted toward the contact, as in case b) of Fig.|6j 
produces the conductance characteristics of the simple 
constriction plus a defect peak at the Dirac energy which 
does not reach the full quantum of conductance. The 
heigth of the peak lowers when shifting the defect more 
distant from the nanoribbon within the graphene lead. 

The fact that the contacts work as diffraction barriers 
decouples the central device from the leads, and hence the 
vacancy state from the leads states. The transmission 
by the vacancy mode through a double barrier can be 
modeled by the standard Lorentzian approximation [1 : 



T{E) 
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r2 [{e-ErY + v^ 



(5) 



where Er is the resonance energy and F = (Fi + F2)/2, 
where Fi and F2 measure the coupling of the resonant 
state to the leads. For a symmetric position of the 
functionalized atom with respect to the left and right 
contacts, like in the positions indicated as a), c) or d) 
in Fig. [6j the double barrier is symmetric. We have 
Fi = F2 = F and the maximum transmission FiF2/F^ 
is one. For an asymmetric position of the functionalized 
atom, like in the position b) in Fig. [6j the double barrier 
is asymmetric and the maximum transmission is smaller 
than one. The maximum transmission in the peak tends 
to zero when the ratio F1/F2 tends to zero or to infinity, 
i.e. when shifting the defect more and more within the 
graphene lead. 

Let us emphasize also that in a metallic nanotube func- 
tionalized by an OH group, a state (for example a tt — tt* 
bond) is removed from a conduction level and the con- 



ductance becomes at the level of the removed state 
[28 . Instead, when functionalizing within diffraction bar- 
riers like in the present case, the conductance is enhanced 
rather than destroyed. 

Let's now discuss experimental implications of our 
findings on quantum dots and suggest new experiments 
as well as possible technological devices. 

We find that the computed characteristics of the geo- 
metric quantum dot (Fig. [5| looks very similar to what 
has been effectively measured in a recent experiment 
done on a graphene quantum dot (see Ref. [IJ). The 
present work gives important insight in the role of diff- 
fraction to produce the electron confinement in the dot, 
in particular for energies close to the Dirac point. 

More generally we note that the control of the char- 
acteristics of existing quantum dot is always a difficult 
task and the graphene quantum dot discussed here offers 
new possibilities. For instance, using appropriate elec- 
tric gates it is possible to adapt the wavelength of 
the incident and outcoming electrons in the half planes. 
This will modulate the characteristics of the barriers that 
couple the quantum dot to the leads. It is then possible 
to go from a small to a large barrier transparency, and 
thus to study a transition from a localized to a metallic 
regime. 

Finally the functionalized graphene constriction could 
have applications in spintronics. Indeed, with a spin po- 
larized defect, the quantum dot allows the production of 
spin polarized currents. Control over the spin of the de- 
fect by a magnetic field could provide control over the 
spin current providing a molecular spin valve. 

To conclude our work establishes an important link 
between nanoelectronics and subwavelength optics, that 
could be experimentally tested thanks to the remar- 
quable properties of graphene. The concept of "elec- 
tronic diffraction barrier" is central to the understanding 
of transport properties of subwavelength graphene nan- 
odevices. It allowed us in particular to get insigth in the 
properties of quantum dots and to propose a new kind 
of technological device for molecular spin valves. This 
shows the high potential of the bridging between sub- 
wavelength optics and nanoelectronics with graphene as 
a constituent material. 
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